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Summary 

The feasibility of using satellite data for climate research over 
the Greenland ice sheet is discussed. In particular, we 
demonstrate the usefulness of Advanced Very High Resolution 
Radiometer (AVHRR) Local Area Coverage (LAC) and 
Global Area Coverage (GAC) data for narrow-band albedo 
retrieval. Our study supports the use of lower resolution 
AVHRR (GAC) data for process studies over most of the 
Greenland ice sheet. Based on LAC data time series analysis, 
we can resolve relative albedo changes on the order of 2-5%. 
In addition, we examine Scanning Multichannel Microwave 
Radiometer (SM MR) and Special Sensor Microwave Imager 
(SSM/I) passive microwave data for snow’ typing and other 
signals of climatological significance. Based on relationships 
between in situ measurements and horizontally polarized 19 
and 37 GHz observations, wet snow regions are identified. 
The wet snow regions increase in aerial percentage from 9% 
of the total ice surface in June to a maximum of 26% in 
August 1990. Furthermore, the relationship between bright- 
ness temperatures and accumulation rates in the northeastern 
part of Greenland is described. We found a consistent 
increase in accumulation rate for the northeastern part of the 
ice sheet from 1981 to 1986. 

1. Introduction 

1.1 Greenland Ice Sheet - Overview 

The Greenland ice sheet is one of the major ice 
sheets in the world. It is the largest one in the 
Northern Hemisphere and is generally thought to 
have survived throughout the Quaternary. The 
ice sheet presently plays an important role in 


the hemispheric circulation through its topography 
and strong temperature gradient towards the 
northwestern Atlantic. The present ice sheet has 
a surface area of 1.75- 10 6 km 2 and a volume of 
2.65- 10 6 km 3 . These correspond to 11 and 8% of 
the global glacier surface area and the ice volume, 
respectively (Ohmura et al., 1991). The ice volume 
is also equivalent to 6.7 m of water relative to the 
present global sea level. The ice sheet has a typical 
shield shape climbing relatively steeply from sea- 
level to 1000 m and then gently to the highest 
point of 3210 m a.s.l. at 72.3° N and 38.0° W, 
where the thickness of ice is estimated at about 
3000 m. The average thickness of the ice sheet is 
1515m. 

The Greenland ice sheet is a remnant of the 
Pleistocene ice sheet. During the last ice age the 
volume was estimated at 2.92 to 5.59* 10 6 km 3 
(Ohmura et al., 1991), suggesting that the surface 
altitude was higher in the past. Presently the 
equilibrium line of the ice sheet lies at 1 800 m a.s.l. 
at the southern end, and it descends to 750 m a.s.l. 
in the northern Greenland (Ambach and Kuhn, 
1985). The accumulation and ablation area occupy 
82 and 18%, respectively. Within the accumulation 
area the dry snow line is found at 3100m in the 
south descending to 1650 m in the north, with the 
surface area of the dry snow zone being approxi- 
mately 50% of the entire ice sheet. 
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1.2 Greenland Ice Sheet - Climatology 

Snow and ice are key variables in the global 
climate system. They can influence the global heat 
budget through their regional feedback mechanism 
in the varying exchange of heat, moisture and 
momentum between the surface and the atmosphere 
(Dickinson et al., 1987). Therefore, accurate infor- 
mation on cryospheric variations (e.g. large scale 
surface albedo increases or decreases due to change 
in dry snow area extent) are essential for the 
prediction of future climate change. Information 
on the response of the ice sheet to climate warming 
is of crucial importance for reliable projection of 
future sea level. As the Greenland ice sheet is 
located in a relatively warm climate, it is particu- 
larly vulnerable to climate change. Present com- 
puter model simulations of the greenhouse scenario 
predict a regional cooling for the Greenland ice 
sheet, contrary to the general warming trend 
over the Arctic sea ice areas (Manabe et al., 1992). 
Trends for the annual temperature computed 
from land surface stations for the period 1961 to 
1990 show a slight cooling for the southern half 
of Greenland (Chapman and Walsh, 1993). How- 
ever, our present knowledge of the mass and 
energy exchange, and its sensitivity to environ- 
mental conditions, is insufficient to make reliable 
quantitative estimates of future changes associ- 
ated with the greenhouse warming. An observed 
0.23 m/year thickening of the Greenland ice sheet 
south of 72° N was reported by Zwally (1989) 
based on satellite radar altimeter measurements. 
The northern part of the ice sheet was not covered 
by satellite altimeter measurements and at this 
point we don’t know if the Greenland ice sheet 
north of 12° N is growing or shrinking. 

Annual total precipitation and annual accumu- 
lation on the Greenland ice sheet were evaluated 
by Ohmura and Reeh (1991). The analysis is based 
on accumulation measurements of 251 pits and 
cores obtained from the upper accumulation zone 
and precipitation measurements made at 35 
meteorological stations in the coastal region. The 
mean annual precipitation for all of Greenland 
was found to be 340 mm water equivalent, with an 
increasing gradient from north to south. The 
amount of precipitation is regulated primarily by 
atmospheric conditions, such as stability, water- 
vapor content and circulation. The winter 
circulation is strongly dominated by two semi- 


permanent cyclones, the Baffin Bay low to the 
west and the larger Icelandic low to the southeast. 
The ice sheet is located under the weak saddle 
between the two depressions. The summer circu- 
lation is dominated by the pressure ridge extend- 
ing from the northeast towards the center of 
the ice sheet. Changes in the surface radiation 
balance and/or in the cloud amount over the ice 
sheet due to global climate change may alter the 
overall pressure distribution and consequently 
the precipitation. It is beyond the scope of this 
paper to derive precipitation amounts, but the 
aerial extent of dry snow area and its seasonal and 
annual variations as derived from satellite data 
will help to understand the present situation. 

1.3 Greenland Ice Sheet - Research Camp 

Scientists from the Swiss Federal Institute of 
Technology (ETH) in Zurich established a perma- 
nent research camp at the equilibrium line of the 
Greenland ice sheet near Jakobshavn (69° 34’ N, 
49° 17’ W) in spring of 1990 (Fig. 1). Climatological 
measurements were carried out to study the energy 
exchange and the mass balance above, at, and 
below the ice surface during April-September 
of 1990 and 1991. Furthermore, the structure of 
the atmospheric boundary layer within the kata- 
batic wind regime on the ice sheet was recorded 
with daily radiosonde profile measurements. Pre- 
liminary results from the two field seasons are 
summarized in progress reports (Ohmura et al., 
1991, 1992). The research station has been trans- 
ferred in summer 1992 to the University of Colo- 
rado at Boulder to continue the field experiments 
for three additional seasons. 

Micro-climatological measurements at a single 
point are necessary and essential for the study of 
the magnitude of the different energy fluxes such 
as the turbulent fluxes and radiation balance at 
the ice sheet surface. However, a single point 
measurement is not sufficient to discuss the clima- 
tology of the entire Greenland ice sheet. Therefore, 
we propose the use of multispectral satellite data 
to derive certain surface parameters which are 
relevant to the surface climatology. The purpose 
of this paper is to demonstrate the feasibility of 
satellite remote sensing as a technique for climate 
studies over large homogeneous snow and ice 
covered areas. The ultimate goal of the project is 
the statistical analysis of the regional and seasonal 
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Fig. 1. Overview of the Greenland ice sheet. 
The outline of the ice sheet is shown with a 
dotted line, and the land is depicted with a 
solid line. The location of the joint research 
camp of the Swiss Federal Institute of Tech- 
nology and the University of Colorado at 
Boulder is shown at the west side of the ice 
sheet near the coastal village of Jakobshavn. 
The labels A through K give the locations 
of the SMMR brightness temperature time 
series as given in Fig. 14 and 15 


variations of surface parameters. This statistical 
analysis over a period of ten years is needed to 
improve the present knowledge and understanding 
of the Greenland ice sheet surface climatology. 

2. Method 

2.1 Multispectral Satellite Data 

Past and ongoing research have shown that several 
polar surface properties can be mapped through 
a combination of multispectral satellite and aircraft 
data (Steffen et al., 1993; Haediger et al., 1993; 
Williams et al., 1991). These properties include: 
spectrally integrated surface albedo, surface tem- 
perature, shortwave and longwave radiation 
balance, and snow properties. All of these param- 
eters can be linked to the surface energy balance. 
The Greenland ice sheet, due to its large size and 


homogeneous surface, is an ideal test area for the 
development of remote sensing techniques in 
climate research. Unlike sea ice and land areas, 
the three major surface types of the ice sheet 
(glacier ice, wet snow, and dry snow) extend over 
several tens of kilometers at a minimum, which 
makes the application of low resolution satellite 
data feasible for surface type classification. 

The types of satellite imagery used for this study 
cover a range of spatial resolutions and wave- 
lengths. For the visible and thermal infrared 
wavelengths we have used Advanced Very High 
Resolution Radiometer (AVHRR) Local Area 
Coverage (LAC: 1.1 km at swath nadir) and Global 
Area Coverage (GAC: 4 km at swath nadir) data. 
For the passive microwave region, Nimbus-7 
Scanning multifrequency microwave radiometer 
(SMMR) data (10/25/80-8/20/87), and DMSP 
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Fig. 2. Processing of multi-spectral satellite data in combination with in situ measurements for the retrieval of the surface 
parameters relevant to climate process studies 


special sensor microwave imager (SSM/I) data 
(7/9/87-9/30/90) were used. They are presently 
available for both polar regions, including the 
Greenland ice sheet, on CD-ROM from the 
National Snow and Ice Data Center (NSIDC), 
Boulder, Colorado. 

In order to retrieve surface parameters relevant 
for climate process studies from satellite imagery, 
a number of sensitivity analyses and intercom- 
parisons with ground measurements are needed. 
These are described below (Fig. 2). 

22 Processing of AVHRR Data 

The AVHRR data were received as raw data 
that have been quality controlled and assembled 


into discrete data sets (level lb data) at the 
Naval Oceanic Atmospheric Research Laboratory 
(NOARL). AVHRR LAC and GAC channel 1 
(0.58 /im-0.68^m) and 2 (0.725 //m-1. 10 gm) 
images have been calibrated over the entire Green- 
land ice sheet using pre-launch calibration coeffi- 
cients to obtain percent albedo values. The cali- 
bration technique follows the description given in 
the NOAA users guide (Kidwell, 1991) and in- 
cludes the non-linear corrections. In order to 
correct for the geometric distortions caused by the 
curvature of the earth and the scanning geometry 
the images have been geo-registered and mapped 
to a polar stereographic projection using a navi- 
gation code developed by Baldwin et al. (1993). 
The accuracy of the navigation routine is at most 
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6 to 8 km due to errors in the time tags. The 
images have been fitted to a map grid to reduce 
this error to less than 1 pixel. The CPU times 
necessary to calibrate and navigate the LAC and 
GAC data for the entire Greenland ice sheet are 
5 min for a GAC image (325 x 625 pixels) and 84 
min for a LAC image (1500 x 2500 pixels) on a 
IBM 6000 workstation. On a SUN II workstation 
13 min for GAC and 188 min for LAC are needed 
to process the satellite imagery. 

Our present computing facility does not permit 
the processing of a large number of LAC scenes 
due to extensive computing time necessary to 
calibrate and navigate the AVHRR LAC images 
for the entire Greenland ice sheet. Thus, it is 
desirable to use the GAC images for long time 
series analysis. In section 3 we present a pilot 
program to test the feasibility of using AVHRR 
GAC images for long term climate studies over 
the Greenland ice sheet. 

2.3 Processing of SMMR and SSM/I Data 

The SMMR is a ten-channel, five frequency linearly 
polarized passive microwave radiometer system. 
The instrument measures surface brightness 
temperatures at 6.6, 10.7, 18.0, 21.0 and 37.0 GHz. 
Vertical and horizontal polarization are provided 
for each frequency. The SSM/I operates at four 
frequencies, namely 19.35, 22.24, 37.0, and 85.0 GHz. 
Both polarizations are provided for each frequency 
except 22.24 GHz. For this study we used the 19 
and 37 GHz frequencies which have effective field 
of view dimensions of 69 x 43 km and 37 x 28 km, 
respectively (along track x cross track). The 
SMMR and SSM/I data provided by NSIDC are 
gridded in polar stereographic projection in 25 x 
25 km grid cells. 

The passive microwave data are read from the 
CD-ROM, and all points within the rectangle 
whose edges contain the outermost limits of the 
Greenland boundary are extracted. Then the 
ocean mask (which is provided on each CD- 
ROM) is applied. This procedure is followed for 
all of the available dates, and it yields a set of 
images 60 x 109 pixels in size for each channel 
and each day of SSM/I and SMMR coverage. A 
total of 1611 images for each SMMR channel and 
1186 images for each SSM/I channel have been 
extracted. 

Some applications, such as snow mapping 
(described in section 3.2) or SMMR/SSM/I cor- 


relations which are described below, require that 
the data be left in image format. For such applica- 
tions an ice mask was developed from a vector file 
of the ice sheet boundaries which was digitized 
from the Geological Survey of Greenland (GGU) 
geological map (scale = 1:250,000). This mask is 
applied to the images and the result was a set of 
images of the Greenland ice sheet with no land- 
contaminated pixels. 

Other applications such as time series investi- 
gations require that data points from the same 
location(s) be extracted from all of the images. 
The first such application is the retrieval of time 
series of brightness temperatures along a horizontal 
transect across Greenland. The transect chosen in 
the first stage of the data analysis is the one that 
passes through the ETH/CU research station 
(Fig. 1). In addition to providing information 
about the brightness temperatures along a west/ 
east profile across the ice sheet, this time series 
provides a means of identifying days of bad data. 
The brightness temperatures between 318-325° E 
longitude are relatively stable with little fluctua- 
tion. Consequently, days that show very high 
deviation from the normal (greater than 4 sigma) 
in this longitude range are flagged as "bad” days, 
and are linearly interpolated from the nearest 
"good” days. 

From the transect, the pixel that corresponds to 
the station location is extracted, and a time series 
of data specifically for the station is produced. In 
addition, the horizontally polarized gradient ratio 
is calculated using the following formula: 

GR = (19 H - 37 H)/( 1 9 H + 37 H) (1) 

where 19 and 37 correspond to the channel 
frequencies in GHz, and H indicates horizontal 
polarization. Similar time series can be produced 
for any location on the ice sheet, but the climate 
station is chosen initially so that comparisons to 
in situ data can be made. 

One final item in the processing of the passive 
microwave data pertains to the relationships 
between SMMR and SSM/I data. An extended 
time series is important for the assessment of 
climatological trends, but in order to effectively 
combine the SMMR and SSM/I data sets, a 
comprehensive understanding of how the two are 
related must be obtained. Figure 3 shows the 
relationship between the 4 common channels on 
the two instruments for the overlapping days of 
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Fig. 3. SMMR and SSM/I brightness temperatures for the overlapping days in 1987 (July 9-August 20) at 72° N and 37° W 
longitude 



SSMI 


Fig. 4. Relationship between SMMR 
and SSM/I data over the entire ice 
sheet for the full period of overlapping 
coverage 


operation in the highest area of the ice sheet (72° N 
and 37° W). A large discrepancy (nearly 10 degrees 
consistently) is apparent between the 19 GHz 
SSM/I channels and the 18 GHz SMMR channels. 
Similar discrepancies, though not as pronounced, 
exist for the 37 GHz channels. 


These relationships were addressed for Antarctica 
by Jezek et al. (1991), but they have not yet been 
defined for the Greenland ice sheet. For all of the 
overlapping days, July 9, 1987 through August 20, 
1987, SSM/I 19 GHz vertically polarized bright- 
ness temperature measurements are compared to 
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the corresponding SMMR observations (18 GHz 
vertically polarized), and the relationships are 
assessed. The results, which correspond well to 
Jezek et al. (1991) results are shown in Fig. 4. The 
small discrepancy can most likely be credited to 
the fact that Antarctica has much colder regions 
than Greenland. Consequently, the data points 
span a much broader range of temperatures, so 
the behavior of the correlation curve is better 
defined at the low end. Furthermore there is 
simply more data from Antarctica because of its 
size. The similarity between the two regression 
lines, however is encouraging. 

Based on this comparison, the SMMR time 
series can be extended to include the SSM/I data, 
for a data set that currently spans over 13 years. 
Similar relationships need to be determined for 
the other channels. 

3. Results 

3.1 Narrow-Band Albedo of Snow 

The variation of solar radiation absorbed and 
reflected by the earth is a key factor in the 
understanding of climate change. If the solar 
radiation reflected by the earth were to increase, 
a trend toward the extension of the winter season 
would result, as can be shown in a simple energy 
feedback mechanism (Budyko, 1974). The snow 




Longitude (degrees west) 


Fig. 6. West-east transect plots of narrow-band planetary 
albedo (0.58-0.68 /im) for AVHRR Local Area Coverage 
(LAC: dotted curve) and Global Area Coverage (GAC: line 
curve) data at 69.6° N of the Greenland ice sheet. Four 
regions for June 06, 1990 are represented. Region A shows 
the bare ice, region B the wet snow, region C the dry snow, 
and region D the clouds 



Fig. 5. Hemispheric spectral albedo for very new 
snow (a few hours old), new snow (12 days old), 
onset of melt, and melting snow surface in the spectral 
range 300-2500 nm measured at the ETH/CU re- 
search camp in 1991 
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cover is one of the most reflective naturally 
occurring materials, and therefore, an under- 
standing of its visible solar reflectance is of great 
importance. The spectral albedo of different snow 
types at the ETH/CU research camp was measured 
for the purpose of calibrating and interpreting 
satellite derived values (Fig. 5). All measurements 
were carried out at solar zenith angles between 
50° and 68°. The snow grain diameter varied for 
the different measurements: very new snow 0.5 mm; 
new snow 0.1 -0.2 mm; onset of melt and wet 
snow 1-5 mm. 

AVHRR LAC and GAC data were processed as 
described in section 2.2 to study the differences in 
pixel resolution for albedo retrieval over the ice 
sheet. Figure 6 shows transect plots of the planetary 
albedo values for both the LAC and GAC channel 1 
(0.58-0.68 /an) across the ice sheet at 69.6° N. 
Four different regions can be distinguished: 

Region A: Bare Ice and Slush 

In this region the majority of the snow has melted 
and bare ice remains. Surface lakes tend to form 
on top of the bare ice, and some slush also exists. 
The agreement between the LAC and GAC images 
is low in comparison to the dry snow regions. 
Such a result is expected because the melt lakes 
are within the 1 km foot-print of the LAC data but 
the GAC footprint contains information from the 
surrounding area as well (Table 1). 

Region B: Wet Snow 

In the wet snow region, the agreement between 
the LAC and GAC data is much improved. 
Although, small scale variation is not detected 
with the GAC data, the general gradual increase 
of albedo is well represented (Table 1). 


Table 1. Statistical Results of the Comparison Between LAC 
and GAC Data from June 06, 1990. Values are given for a 
region 300km 2 region 


Region 

Mean 

difference % 

Std dev of 
difference 

Correlation 

Bare ice 

-2.05 

3.06 

0.890 

Wet snow 

-0.16 

0.89 

0.970 

Dry snow 

-0.04 

0.26 

0.998 

Clouds 

0.06 

2.21 

0.997 


Region C: Dry Snow 

Overdry snow where the albedo remains relatively 
constant, the agreement between the LAC and the 
GAC images are very high, with mean differences 
practically zero. In this region the albedo is very 
homogeneous and thus very little difference between 
the two resolutions is expected (Table 1). 

Region D: Clouds 

Detection of clouds tends to be difficult since 
cloud reflectance is similar, if not identical, to 
snow reflectance. However, thick overcast regions 
show noticeable fluctuations in reflectance (5%). 
In these regions the agreement between the LAC 
and the GAC images decreases slightly (Table 1), 
but due to the large size of the cloud features (4 km 
or larger) the GAC data is able to detect the slight 
variation in reflectance due to the presence of 
clouds as well as the larger decrease in reflectance 
from the cloud shadows. The ability of the GAC 
data to detect the cloud shadows could be a useful 
tool in cloud discrimination over the Greenland 
ice sheet. 

Thus, the same information for wet snow, dry 
snow, and clouds can be obtained by using either 
high or low resolution data. Using the GAC data 
during summer months, when the melt region is 
relatively large, should be adequate. When. the 
melt season is just beginning, the lower resolution) 
data will not be able to accurately detect the 
melting since the pixel will be contaminated by 
dry snow. Higher resolution data such as LAC 
will be needed at the margin of the ice sheet in 
such cases. In addition, to study the surface lakes 
that form during the melting season, high resolution 
data with pixel size less than 2 km (average lake 
diameter) will also be needed. 

To obtain snow surface albedo maps it is 
necessary to correct the AVHRR satellite data for 
the intervening atmosphere. Such a correction can 
be accomplished using a linear relationship between 
clear sky planetary albedo and surface albedo 
(Koepke, 1989), or by using a radiative transfer 
model, such as LOWTRAN or 5S, to correct for 
scattering and absorption. Although presently 
there is no consensus among scientists as to which 
atmospheric radiative transfer code is best, most 
codes tend to disagree only for large optical 
depths of aerosols and large off-nadir angles (60° 
or more) (Royer et al., 1988). The 5S code developed 
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by Tanreet al. (1990) provides reasonably accurate 
atmospheric modeling and surface reflectance 
retrieval (Teillet, 1992). The overall scheme involves 
correcting the narrow-band planetary albedo as 
viewed by the sensor to the narrow-band surface 
albedo, taking into account illumination, sensor 
view angle, and atmospheric effects. 

To run the code, illumination and observation 
geometries for each pixel are provided as input. 
The water vapor input can be calculated from 
the AVHRR channels 4 (10.30-1 1.30 /im) and 5 
(1 1.50- 12.50 //m) or it can be provided by radio- 
sonde data that has been collected at the ETH/ 
CU camp during the summers of 1990 and 1991 
(Haefliger et al., 1993). Radiosonde data are also 
available from the summit of the Greenland 
ice sheet and from the coastal station Gothavn 
and will serve to derive an “average atmosphere” 
for the data available in 1990 and 1991. Further 
parameterization is needed to extrapolate for 
other years. The ozone concentration is taken to 
be that of subarctic summer. Aerosol content can 
vary significantly in Arctic regions (Lindsay and 
Rothrock, 1993) and can result in large errors in 
deriving the surface reflectance. Unfortunately, 
optical depths of aerosols over the Greenland ice 
sheet are not known and standard atmosphere 
assumptions must be made. Using the 5S code, it 
is possible to retrieve surface reflectance by the 
following expression: 

Pi = 100 {Ai + B,) ( 2 ) 

(\00 + (A i + Bi)Si) 

A . = — t B j — HzJOO^gf ?) (3) 
V s T„COS0 s V. 

where the index j stands for the AVHRR channel 

(i = 1,2). 

Pi = surface reflectance for channel 1 or 2 in 
percent 

p° = reflectance at the sensor 
Patm = atmospheric reflectance 
S = atmospheric spherical albedo 
t g = total gas transmittance 
t s = total scattering transmittance in solar 
direction 

t v = total scattering transmittance in sensor 
direction 

d s = solar distance 
6 S = solar zenith angle. 


Table 2. Input Parameters for the 5S Code Run for the Location 
of the ETH/CU Camp, May 23, 1991 


Solar zenith angle 

49. 0 C 

Sensor zenith angle 

29.40 

Solar azimuth angle 

285.57 

Sensor azimuth angle 

243.69 

Atmospheric profile 

(i) H 2 0 = 0.62 gem 2 

(ii) 0 3 = subarctic summer 

Aerosol model 

Continental 

Aerosol level 

Visibility = 8 km 

Spectral bands 

NOAA-ll AVHRR channel 1 

Surface reflectance 

Dry snow 


These quantities, which are necessary to calculate 
the surface reflectance, are obtained by running 
the 5S code. 

A three day time series for the AVHRR LAC 
channel 1 was atmospherically corrected using the 
in situ radiosonde profile data (See Table 2 for 5S 
inputs). Figure 7 illustrates the temporal changes 
of narrow-band surface albedo for the dates May 
23, June 01, and June 08, 1991, along a transect 
at 69.6° N (see also Fig. 1). This time series shows 
a continual decrease in albedo over time. On May 
23, 1991, the snow cover was dry and the land was 
also covered by snow with many rocky areas 
showing through. By June 01, melting had occurred 
and the surface reflectance dropped considerably, 
and by June 08 the reflectance has decreased 
further by 5%. The ability of the sensor to detect 
melt enables us to determine the aerial extent of 
the transition regions along the perimeter of the 
ice sheet. Figure 8 shows the AVHRR channel 1 
imagery for the Jakobshavn area on June 8, 1991 
(see also Fig. 1 for location). Cloud shadows 
running north-south on the left side of the image 
(snow covered ice sheet) can be identified. The 
snow free areas of the ice sheet along the margin 
(bare ice) stand out against the higher reflectance 
of wet snow areas. 

The narrow-band surface albedo as derived 
from AVHRR channel 1 and corrected with the 5S 
radiative transfer model for atmospheric scatter- 
ing, and absorption was compared with in situ 
measurements at the ETH/CU station. The date 
used for this intercomparison is 23 May 1991 
under clear sky conditions. The in situ hemispheric 
spectral albedo, integrated over the wavelength 
region 580-680 nm, is 91.6%, and the narrow- 
band albedo derived from AVHRR channel 1 is 
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Fig. 7. Narrow-band surface albedo 
(0.58-0.68 fxm) time series showing 
the melt region at 69.6 N. The time 
series is from May 23, June 01 and 
June 08, 1991 



Fig. 8. AVHRR channel 1 imagery for 
the Jakobshvan area of the Greenland 
Y i ce sheet, June 8, 1991. Location and size 

of this imagery is given in Fig. 1. Bare 
ice, wet snow, dry snow and cloud shadows 
A can be identified by different gray levels 

dry snow | cloud shadows (reflectance) 


90.2%. A recent intercomparison with the same 
input values but using the LOWTRAN for the 
atmospheric correction gave a narrow-band 
albedo value of 93.3% (Haefliger et al., 1993). 
This preliminary analysis shows that with in situ 
atmospheric temperature and humidity profile 
data as input variables for the radiative transfer 
model, narrow-band surface albedo can be derived 
within one and two percent accuracy from AVHRR 


channel 1. Furthermore, the comparison suggests 
that the standard aerosol and ozone concentrations 
give reasonable results for this case study. However, 
comparisons for different seasons and a larger 
number of cases are needed to derive statistically 
significant values. 

It should be noted that altitude dependence was 
not considered in these computations, as we were 
only interested in relative changes of the albedo. 
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In order to derive narrow-band albedo maps of 
the entire Greenland ice sheet, this correction is 
essential. However, the discussion of albedo maps 
for the entire Greenland ice sheet is beyond the 
scope of this paper. 

3.2 Dry Snow/Wet Snow Classification 

The passive microwave brightness temperature 
T b is, according to Rayleigh Jeans approximation, 
a product of the emissivity e and the mean 
physical temperature Tof the snow layer from 
which the radiation is emanating. 

T„ = sT. (4 ) 

Radiative transfer model calculations show that 
the penetration depth for dry snow depends 
mainly on snow density, the radius of the snow 
grains particles, and the physical temperature of 
the snow medium (Zwally, 1 977). The penetration 
depth of passive microwave radiation is defined 
as the inverse of the extinction coefficient. It is a 
function of frequency, grain size, physical temper- 
ature of the snow layer, and density. According 
to Ulaby et al. (1986) the penetration depth (d) 
for snow can be approximated by: 

d = (A/2n)(ftf/7'), (5) 

where X is the wavelengths of the radiation in 
meters, e' is the relative permittivity which is a 
constant for snow (e' = 3. 1 5), and e" is the dielectric 
loss factor. The penetration depth was calculated 
according to Eq. 5, and the values are given in 
Table 3 for two frequencies and two different 
temperatures. 

For melting snow, the penetration depth is 
strongly reduced by increasing liquid-water content 
in the snow pack. For wet snow the volume 
scattering is absent due to strongly reduced pene- 
tration, which leads to the near-black body emis- 
sion, especially at vertical polarization (Matzler 


and Hueppi, 1989). Thus it can be seen that the 
microwave emission emanates from different depths 
for different snow types (dry and wet snow). At 
frequencies above 20 GHz, scattering becomes the 
dominant component of the total extinction loss 
of the medium. This means that the emissivity of 
the snow depends strongly on the snow grain size, 
and increases with decreasing grain radius (Sri- 
vastav and Singh, 1991). Based on these theoretical 
considerations we would expect a large variation 
of the brightness temperature over the Greenland 
ice sheet which should yield information on the 
different snow facies associated with melt and 
accumulation rates. 

Figure 9 shows a brightness temperature time- 
sequence profile across Greenland from 310°- 
335° E, at 69.6° N for the SSM/I 19 GHz vertical 
polarization. The record is for January 13, 1988 
through March 31, 1991 at a daily interval. The 
three main melting periods at the margins of the 
Greenland ice sheet are clearly shown in the 
graph. The brightness temperature increase is due 
to near black-body emissions of the wet snow 
cover. For the higher regions of the ice sheet 
between 320° and 325° E, where the air temperature 
is always below freezing (dry snow regions), the 
T b is quite constant. Figure 10 shows two selected 
brightness temperature cross sections to demon- 
strate the difference between the summer-melt 
condition and the winter condition. The profile 
consists of 33 SSM/I pixels (25 x 25 km). The 
increase of the winter T b at the east coast of 
Greenland is due to the close proximity of land 
(spill-over) and can be explained by large 19 GHz 
field of view (69 x 43 km) binned to a 25 x 25 km 
grid. 

Equation 4 shows that T b is related to the 
physical temperature, which makes snow type 
classification based on a single channel impossible 
for large regions extending over different tem- 
perature regimes. The gradient ratio (GR) between 


Table 3. Microwave Penetration Depths for Dry Snow with a Grain Size Radius 
Mean Snow Temperatures of - 1 C and - 20 C 


of r = 0.5 mm at 19 GHz and 37 GHz 


for a 


Frequency 

Dielectric loss 
factor at — 1 °C 

Penetration 
depth (m) 

Dielectric loss 
at — 20 C 

Penetration 
depth (m) 

19GHz 
37 GHz 

6.8 10 3 
2 10 2 

0.66 

0.10 

1.8 • 10 _ 3 
2.7- 1 0 _ 3 

2.5 

0.83 
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Variation of Brightness Temperature at 69.6 N lat (19v) 



Fig. 9. Passive microwave brightness 
temperature time series across Green- 
land for SSM/1 19 GHz vertical 
polarization. The time series runs 
from January 1, 1988 through March 
31, 1991. The dominant three peaks 
along the ice margins demonstrate 
the T b increase due to wet snow 



Fig. 10. Same as Fig. 9 for two selected days in summer 
and winter 


19 and 37 GHz (Eq. 1) shown in Fig. 11 for three 
consecutive years at the ETH/CU station reduces 
the snow temperature dependence. These values 
were compared with automatic weather station 
data collected within 10 km of the station by the 
Geological Survey of Greenland (unpublished 
data, personal communication Henrik Thomsen, 
GGU). For all three years the sharp GR increases 
in spring (e.g. April 4, 1988) coincide in time with 


the air temperatures increases above the zero 
degree centigrade on the very same day. The sharp 
decreases of the GR value at the end of the 
summer (e.g. September 15, 1988) coincide also 
with sub-zero air temperature readings from the 
automatic weather station. The air temperature 
measurements from the ETH/CU research station 
and the corresponding GR values as derived from 
the SSM/I satellite data are shown in Fig. 12. It is 
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Jan-88 Jul-88 Jan-89 JuI-89 Jan-90 Jul-90 Jan-91 Jul-91 



Jan-88 Jul-88 Jan-89 Jul-89 Jan-90 Jul-90 Jan-91 Jul-91 


Fig. 1 1 . SSM/I gradient ratio, ! 9 GHz and 37 GHz 
horizontal polarization, at the location of the ETH/ 
CU research station on the Greenland ice sheet 
from January 1988 through July 1991 


May-90 Jun-90 Jul-90 Aug-90 Sep- 90 Oct-90 Nov-90 Dec-90 



May-90 Jun-90 Jul-90 Aug-90 Sep- 90 Oct-90 Nov-90 Dec-90 


Fig. 12. SSM/I gradient ratio, 19GHzand 37GHz 
horizontal polarization, and air temperature at 
the location of the ETH/CU research station for 
summer 1990 


rather surprising how well the air temperature 
correlates with the satellite measurements (25 x 
25 km gridded data). The correlation is negative 
due to the different penetration depths of 19 and 
37 GHz. An increase in surface temperature effects 
the 37 GHz frequency more than the 19 GHz 
frequency and consequently the GR value decreases. 


The same response, but with negative sign, can be 
seen with a surface cooling. 

The GR value seems to be a good classifier for 
the distinction between dry and wet snow. Since 
large GR fluctuations, caused by the temperature 
gradient change within the snow cover, are large 
we have averaged the GR value over 10 days for 
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April 1990 
Wet Snow Areas 


May 1990 
Wet Snow Areas 


June 1990 
Wet Snow Areas 



July 1990 
Wet Snow Areas 




September 1990 
Wet Snow Areas 




Fig. 13. Time sequence of wet snow areas as classified with gradient ratio threshold (19 and 37GHz/H) for each mid month 
period of April through September 1990. Each pixel ( + ) represent an area of 25 x 25 km 


each month and used a threshold of GR = + 0.025 
to distinguish between dry and wet surface snow. 
In Fig. 13 we show the wet snow area of the entire 
Greenland ice sheet for each mid-month period of 
April through September, 1990. In April only a 
few wet snow signatures are classified along the 
west coast of Greenland. A very strong signal is 
observed at 80° N in the region of the Humboldt 
Glacier (see Fig. 1 for location). Melting at such 
high latitudes is questionable that early in the 


season; however, the close proximity of the North 
Water polynya, with air temperatures up to 8 °C 
above normal, could explain this early melting 
(Steffen, 1985). This melt signal is consistent 
throughout the following 4 months. In May the 
melt areas have extended along the margin of the 
ice sheet south of 70° N. The southeast coast of 
the ice sheet shows little melting which could be 
an artifact of the large pixel footprint and the 
steeply rising ice sheet in those regions. In June 
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and July the wet snow areas increase from 9% to 
22% of the total ice area and reach the maximum 
extent of 26% in mid August. The gradient threshold 
method classifies wet snow areas, whereas glacier 
ice has a GR value close to zero. This explains 
the “unclassified” areas along the west coast in the 
vicinity of the ETH/CU research camp for the 
August period. Local observations from helicopter 
flights in that region during August 1990 confirm 
this hypothesis (Ohmura et al„ 1992). In September 
the wet snow areas are reduced to the region of 
the low ice margin as the new snow precipitation 
masks the melt signal. 


3.3 Passive Microwave Brightness 
Temperature Variations 

Additional time series of interest are those at 
different locations near the center of the ice sheet 
where the altitude is high, the snow is dry, and the 
region is the most homogeneous. Some interesting 
trends in brightness temperatures are found in 
some of these areas, and they are discussed below. 
At these locations, instead of extracting the ob- 
servations by means of a transect and then a single 
pixel value (as in the case of the climate station), a 
2x2 pixel region is filtered out of each image and 
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Fig. 14. Time series of brightness temperatures for the full period of SMMR coverage at different locations along the central 
ridge of the Greenland ice sheet (see Fig. 1, A-F) 
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Fig. 15. Time series of brightness temperatures for the full period of SMMR coverage at different locations along the central 
ridge of the Greenland ice sheet, off of the central ridge (see Fig. 1, G-J) 


is spatially averaged. In this way, localized anom- 
alies are smoothed out. 

The locations of the selected areas are shown 
in Fig. 1, and the corresponding time series of 
brightness temperatures from the SMMR 18 GHz 
vertically polarized channel are shown in Figs. 14 
and 15. The SMMR data set was selected over 
SSM/I because it spans a greater period of time. 
Ultimately, the two will be combined. Locations 
A through F are chosen because they span the 
highest sections of the ice sheet for the given 
latitude. Effects due to melt are thus eliminated. 
Additional areas of the ice sheet were selected for 
broader comparisons, but not all are shown in the 
figures. 

A trend is clearly evident in most of the plots 
between the winters of 1 98 1 and 1986, with winter 
brightness temperature increases ranging from 
1 K to nearly 20 K over the 5 year period. We 
believe the trend holds some information about 
climatological features for the various locations 
sampled. It may be due to changes in the physical 
temperature, variations in the snow accumulation 
rate, differences in snow grain sizes over those 
years, changes in snow depth, or a combination 


of variations in several of these features. Each can 
affect the brightness temperature. 

Grain Size 

The relationship between grain size and brightness 
temperature has been clearly demonstrated 
(Srivastav and Singh, 1991). As grain size increases, 
the T b also decreases. If there is snowfall in the 
region, the grain size depends on the atmospheric 
temperature and humidity, but should remain 
relatively constant due to wind action. The grain 
size of the snow underneath will increase due to 
snow metamorphosis (Seligman, 1980). If there is 
no accumulation, then snow metamorphosis will 
also cause grain size increase. In either case, the 
effect on brightness temperature will be a decrease, 
or no change. In fact, the only way grain size of 
new snow can increase the brightness temperature 
is by reducing in size. It is not likely that the 
observed T b increase over several years (Fig. 
14 A) is caused by reduction in grain size at the 
surface of the snow layer. The possibility that 
every year there is a constant decrease in tem- 
perature such that the grains become smaller and 
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the brightness temperature consistently increase 
is minimal. Furthermore, the temperature and 
grain sizes are not correlated over a very broad 
range. For these reasons, the observed trend is not 
attributed to grain size changes. 

Physical Temperature 

It is possible that the observed trend may reflect 
a physical temperature rise since the two are 
directly proportional. However, there is no evidence 
of such a temperature increase from the nearby 
meteorological stations at the north-east coast of 
Greenland (Nord and Danmark Havn). In addition 
there is a large variation in the range of tempera- 
ture increase from one location to another. The 
maximum occurs at 78° N and 32° W, and it is 
1 8 K. At the same latitude, but at 48° W longitude, 
there is a minimal change of 1 K. Furthermore, 
the regions south of 70° N seem to exhibit almost 
no apparent pattern. We would expect, if a rise in 
physical temperature were responsible for the ob- 
served brightness temperature increases, that the 
trends would be much more consistent through- 
out the ice sheet. 

Accumulation Rate 

An increase in accumulation rate results in an 
increase in emissivity (Zwally, 1977), which is 
directly proportional to brightness temperature. 
Therefore, it is possible, that the trends observed 
may be explained by changes in accumulation 
rate. 

As stated in section 3.2, the penetration depth 
for the 19 GHz vertical channel at — 20 °C is 
2.5 m. Winter temperatures averaged from 1951 
through 1960 were found to be — 40 °C in the 
northeastern portion of the ice sheet (Ohmura, 
1987), which is where the trend is strongest. There 
is no published dielectric loss factor for such a low 
temperature, but according to Eq. 5, the penetra- 
tion depth will be even greater than the 2.5 m 
calculated for — 20 °C. 

For locations in which the mean accumulation 
rates are higher than the penetration depth during 
the winter season, snow accumulations do not 
support a sustained brightness temperature 
increase. For extreme cases in these areas, a year 
of very low accumulation might stand out as a 
drop in brightness temperature as the emissivity 
of larger subsurface grain sizes serves to reduce 


the signal. However, a trend of high accumulation 
rates in these areas over several years should only 
result in a rise in brightness temperatures for the 
first one or two years. After that time, the curve 
should level off as the signal contribution of the 
smaller grained new snow increases while the 
contribution of the old snow below diminishes. 

On the other hand, if the mean accumulation 
rate is relatively low in comparison to the pene- 
tration depth, and there is always a significant 
contribution to the signal from old snow below 
the surface, then changes in accumulation rates 
should correspond directly to changes in brightness 
temperatures. High accumulation rates in these 
typically low rate regions will result in effectively 
smaller grain sizes being responsible for most of 
the signal, and will appear as an increase in T b . 
Low accumulation rates will result in effectively 
larger grain sizes due to the strong contribution 
of the old snow below the new, and thus reduce 
the brightness temperatures. In these areas the 
sensitivity of T b varies inversely with the mean 
accumulation rate (Zwally, 1977), in other words 
lower accumulation rate corresponds to higher 
T b and visa versa. 

Mean accumulation rates throughout the ice 
sheet were evaluated by Ohmura and Reeh (1991). 
Their findings indicate the lowest rates, on the 
order of 100-200 mm water equivalent annually, 
are characteristic of the northeast portion of the 
ice sheet. The brightness temperatures in this area 
should be most sensitive to accumulation, and 
they do in fact exhibit the strongest trends (Fig. 

1 5 G). Conversely, the rest of the ice sheet, west of 
the divide, and also south of 70 c N has high 
accumulation rates, and the corresponding bright- 
ness temperature trends are very small. Further- 
more, the brightness temperature trends increased 
with higher latitude (see the temperature changes 
stated in Fig. 14), while coincidentally, the accumu- 
lation rates decreased. 

The interannual variations in snow accumu- 
lation for 78° N and 32° W for 1979 through 1987 
are qualitatively described in Table 4. Similar 
relationships can be determined for any other 
location in the sensitive northeast region. 

Snow Depth 

The argument for snow depth is similar to that for 
accumulation rate. Deep snow on the ice, which 
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Date 


Fig. 16. Combined time series of 
SMMR and SSM/1 data at 78° N 
latitude and 32 c W longitude. The 
offset between the two data sets 
was corrected according to the linear 
fit given in Fig. 4 


Table 4. Relative Accumulation Rates at 78 v N and 32° W 
( NE Greenland ) for 1979 Through 1987 as Determined from 
Fig. 15 


Year 

Comparison to previous year 

1979 

less accumulation 

1980 

less accumulation 

1981 

less accumulation (least of all) 

1982 

more accumulation 

1983 

more accumulation 

1984 

more accumulation 

1985 

more accumulation 

1986 

about the same accumulation 

1987 

same accumulation 


results from high accumulation rates, exhibits a 
higher brightness temperature than shallow snow 
at the same physical temperature. The reason is 
that old subsurface snow from previous years will 
contribute less to the microwave signal if the 
accumulation rate is high. As with accumulation 
rate, the relationship is only significant in regions 
where the depth of new snow is less than the 
penetration depth of passive microwave radiation. 

Figure 16 shows a comprehensive time set of 
18 GHz vertically polarized SMMR data and 
19 GHz vertically polarized SSM/I data. Due to 
overheating, the SSM/I sensor was shut off during 
the winter of 1987, and consequently there is a two 
month gap in the data. The SSM/I data were 


adjusted according to our regression coefficients 
displayed in Fig. 4. As with the SMMR data, 
interannual changes in the brightness temperature 
are apparent in the region with a low snow 
accumulation rate. A decrease in accumulation 
rate during the period of SSM/I coverage is 
indicated by the decrease in brightness tempera- 
ture. 


4. Conclusion 

Little difference between the lower and higher 
resolution AVHRR data is seen for regions of dry 
snow, melt areas, and clouds. Therefore, Global 
Area Coverage (GAC) data can be used for 
long-term snow surface albedo studies over most 
of the Greenland ice sheet. Local Area Coverage 
(LAC) data is found to be necessary during 
transitional periods in the spring, and thus, is 
needed for studies of the onset of melt. 

Small albedo changes of a few percent can be 
detected, which enables us to classify different 
snow types in conjunction with ground based 
spectral measurements. For absolute albedo de- 
termination, in situ radiosonde, ozone, and aerosol 
measurements are required for input into radiative 
transfer models. In addition, a digital elevation 
model is needed to correct for the changes in 
atmospheric path length over the ice sheet. 
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There is a distinct passive microwave signal of 
melting, which can be used to classify wet snow 
areas and identify the onset of melt. This is 
accomplished by using a gradient ratio thres- 
hold of SSM/I 19 GHz and 37 GHz horizontal 
polarization. The seasonal change of wet snow 
extend was analyzed for the months May through 
September, 1990. Wet snow areas are generally 
larger on the west side of the ice sheet due to lower 
topography and consequently higher air tempera- 
tures. The maximum wet snow extent of 26% of 
the total ice area was found in August. The Bare 
ice areas of the margins, where the seasonal snow 
cover has disappeared, are not included in this 
statistics. 

Many features within the snow pack contribute 
to variations in brightness temperature. These can 
be attributed to real physical changes that are 
important for climate sensitivity studies. The issue 
is further complicated by variations of penetration 
depth of the two frequencies used. A better under- 
standing of the physical processes, as well as in 
situ data are needed in order to more accurately 
identify and quantify those features. 

However, the evidence does suggest that in 
some regions of the ice sheet, relative or qualitative 
statements can be made about the rate of accumu- 
lation. This is very encouraging since there is not 
yet an accurate method for quantifying snow 
cover over ice (most of the algorithms to date deal 
with snow over land). The accumulation increase 
as detected for the north-eastern part of Greenland 
during 1981 through 1986 is in general agreement 
with average ice-sheet surface elevation changes 
for the southern part of Greenland as reported by 
Zwally ( 1 989). More work needs to be done before 
such quantitative statements about the snow 
accumulation can be made, but is does appear to 
be feasible. 

While snow typing is easily accomplished with 
AVHRR data during clear sky conditions, it is still 
limited by cloud cover. This limitation can be 
overcome by using passive microwave data, whose 
signal penetrates clouds. However, passive micro- 
wave data is limited by its coarse spatial resolution; 
thus, important sub-pixel phenomena cannot be 
resolved. A combination of the two data sets 
seems promising and will yield a more complete 
description of surface conditions than can be 
obtained by each individually. 
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